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In a recent letter one of us pointed out how differences in preparation procedures for quantum 
experiments can lead to non-trivial differences in the results of the experiment. The difference 
arise from the initial correlations between the system and environment. Therefore, any quantum 
experiment that is prone to the influences from the environment must be prepared carefully. In 
this paper, we study quantum process tomography in light of this. We suggest several experimental 
setups, where preparation of initial state plays a role on the final outcome of the experiment. We 
show that by studying the linearity and the positivity of the resulting maps the experimenter can 
determine the nature of the initial correlations between the system and the environment. 



Quantum information processing promises powerful 
computational methods that surpass the methods of clas- 
sical information processing [TJ [5] . These methods rely 
on taking advantage of quantum parallelism by using 
quantum superposition and quantum entanglement as re- 
sources. In order to implement such a device one must 
have precise control over the system, and isolate it from 
the surrounding environment to preserve coherence. Yet 
realistically, it is nearly impossible to isolate the system 
of interest completely from its surroundings, while having 
a great deal of control. 

With the rising interest in quantum computation and 
quantum information processing, quantum coherence ex- 
periments are performed readily these days, though with 
relatively small systems. One of the major problem with 
these experiment is the loss of coherence due to the inter- 
action between the system of interest and the unknown 
environmental states. The methods for studying the in- 
teraction between the system and the environment are 
given by the quantum theory of open systems. 

The quantum theory of open systems got its start in 
almost fifty years ago with the introduction of dynamical 
maps O 0] due to Sudarshan, Mathews, Rau, and Jor- 
dan. Decades after its conception, the dynamical map 
formalism is finally being tested in the laboratory set- 
ting. The experimental determination of a dynamical 
map is achieved by a procedure called quantum process 
tomography. 

Any experiment, including quantum process tomogra- 
phy experiments, requires a method to prepare the initial 
states of the system at the beginning of the experiment 
[5, 6, 3 13- We study the affects of the preparation pro- 
cedure on quantum systems that interact with an envi- 
ronment. The act of preparation has been neglected from 
the theory of quantum process tomography (and for all 
quantum experiments that interact with a non-trivial en- 
vironment). We investigate this issue for quantum pro- 
cess tomography in detail in this paper. We present sev- 
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eral simple examples to motivate our arguments. 

In previous papers we studied not completely positiv- 
ity of dynamical maps as function of initial correlations 
between the system and the environment [9^ and the role 
of state preparation in quantum process tomography [5] 
(and quantum mechanics in [8|). This paper is an exten- 
sion of those studies. 

In Sec. [T] we review the concepts and the mathematics 
of quantum process tomography and preparation proce- 
dure. In Sec. [n] we construct and discuss several simple 
examples of quantum process tomography for different 
preparations procedures. In Sec. |III[ we follow these ex- 
amples up with analyzing some recent quantum process 
tomography experiment. The example are in the same 
spirit as the ones in the previous papers, therefore al- 
lowing us to compare all of the results. We reproduce 
the examples from previous papers in the Appendix of 
this paper for ease. We have modified the language of 
these examples to fit the language of this paper. And fi- 
nally in Sec. |IV| we discuss how preparation procedures 
differentiate between the outcome of a quantum process 
tomography experiment and its theoretical analogue, dy- 
namical maps, along with our concluding remarks. 



I. A BRIEF REVIEW 
A. Quantum Process Tomography 

Quantum process tomography [TUl E] is the experi- 
mental tool that determines the open evolution of a sys- 
tem that interacts with the surrounding environment. It 
is the tool that allows an experimenter to determine the 
unwanted action of a quantum process on the quantum 
bits going through it. It is an important tool for quantum 
information processing. A state going through a quan- 
tum gate or a quantum channel will experience some in- 
teractions with the surrounding environment. Quantum 
process tomography allows the experimentalist to distin- 
guish the differences between the ideal process and the 
process found experimentally. Therefore it is an impor- 
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tant tool in quantum control design and battling deco- 
herence (loss of polarization). 

The objective of quantum process tomography is to de- 
termine how a quantum process acts on different states of 
the system. In very basic terms, a quantum process con- 
nects different quantum input states to different output 
states: 



input states process output states. 



(1) 



The complete behavior of the quantum process is known 
if the output state for any given input state can be pre- 
dicted. 

The tomography aspect of quantum process tomogra- 
phy is to use a finite number of input states, instead of 
all possible states, to determine the quantum process. 
For instance, to determine a dynamical map, B, we only 
need to know the mapping of the elements of the density 
matrix from an initial time to a final time. 



B 



'(to) Prsit). 



(2) 



The elements of the density matrix linearly span the 
whole state space |7j. 

Experimentally, we do not have the access to the indi- 
vidual elements of the density matrix; we can only pre- 
pare physical states. Thus, a set of physical states that 
linearly span the state space will be sufficient for the ex- 
periment. A state space of dimension d requires (P states 
to span the space. Once the evolution of each these input 
states is known, by linearity the evolution of any input 
state is known (see [1] for detailed discussion). 

Using the set linearly independent input states P^"^\ 
and measuring the corresponding output states Q^"^\ the 
evolution of an arbitrary input state can be determined. 
Let A be the map describing the process, which we call 
process map, and an arbitrary input state be expressed 
(uniquely) as a linear combination ^ ^ p'^'^) p'^"^). The 
action of the map in terms of the matrix elements is as 
follows: 



1131 1141 I15j . direct characterization of quantum dynam- 
ics [Ml [17], selective efficient quantum process tomogra- 
phy^^, and symmetrized characterization of noisy quan- 
tum processes [19". Some of these procedures have been 
experimentally tested ^i2]ji22j|13J|21|251[2a|271[2H]. 



1. Weak coupling assumption 

In every quantum process tomography procedure listed 
above the input states are thought to be pure states. 
There are two advantages of using pure states as inputs. 
First, it is easier to span the space of the system with 
a set of pure state than it is with a set of mixed states. 
Second, pure states are always uncorrelated, which is one 
of the central assumption in every procedure above. This 
is often called the weak coupling assumption |29j . It is 
simply a matter of preparing necessary pure states to 
perform a quantum process tomography experiment. 

What if we depart from this assumption? It is well 
known that that initially correlated states can lead to 
not-completely positive dynamics for the system [301 1311 
I32j . In many recent experiments, the process maps that 
characterize the quantum operations have been plagued 
with negative eigenvalues and occasional non-linear be- 
havior. We now examine the how the two cases are re- 
lated. 

The quantum process tomography procedures we re- 
viewed require the input states be uncorrelated with the 
environment. Just before the experiment begins, in gen- 
eral the state of the system may be correlated with the 
environment. Thus, at the beginning of the experiment it 
is necessary to prepare the initially correlated total state 
into an uncorrelated state. The difference between the 
process maps found from a quantum process tomography 
experiment and the dynamical maps calculated theoret- 
ically is precisely the act of preparation of input states. 
Let us investigate this issue by analyzing the steps in- 
volved in a quantum process tomography experiment. 
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B. Preparation procedure 



With the knowledge of the output states corresponding 
to the input state we find the map by the following ex- 
pression. 



(m) 

^^rr';ss' — / ^ ^:rs ^ r' s' ' 



(3) 



where P*^"^ are the duals of the input states satisfying 
the scalar product 



p(m)t p(ri) _ p(m)* p(n) _ X 

r r ~ ^rs ^rs ~ 



Today there are many variations of the quantum pro- 
cess tomography procedure described above, namely an- 
cilla (entanglement) assisted process tomography jl2| 



In practice, preparation procedures can be very com- 
plicated. Instead of describing many different proce- 
dures, we follow the theory of preparations developed in 
[8]. Below is a brief review of the findings in that paper. 

A preparation procedure is the mapping of a set of 
unknown states into a fix known input state. The most 
general transformation of a quantum state is described 
by a stochastic map [3]. In light of that, a very com- 
plicated preparation procedure due to the aparatus, is 
simply denote by stochastic map, The procedure 

of preparing the mth input state is given as 



P 



S£ 



R 



(4) 



where T is the identity map acting on the state of the 
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environment. In other words, we are assuming that the 
preparation procedure acts only on the system^. 

In general, the state of the system and the environment 
before preparation, p^^ , will be correlated (see [33]). The 
goal of the preparation procedure is to eliminate these 
correlations, as for an ideal experiment the state of the 
system should be uncorrelated with the environment at 
the beginning of the experiment. Therefore, under perfect 
preparation procedure assumption, the right hand side of 
Eq. |4] or the post-preparation procedure state, is of prod- 
uct form. 

In [S] two common preparation methods, that achieve 
product states, were discussed. In this paper we will only 
focus on the method called stochastic preparation. The 
other method, projective preparations will be discussed 
elsewhere. 



Stochastic preparation 



Many quantum experiments begin by initializing the 
system into a specific state. For instance, in the simplest 
case, the system can be prepared to the ground state by 
cooling it to near absolute zero temperature [23 (HI 
[26] . Mathematically, these set of operations axe written 
as a pin map |34j . 



(5) 



where 1 is the identity matrix, acting as the 'trace op- 
erator' {IrsPrs — Tr[p]) and |$) is some fixed state (i.e. 
ground state) of the system. In this procedure, no matter 
what the initial state of the system was, it is "pinned" 
to the final state 

The action of the pin map on a bipartite state of the 
system and the environment. 



e 



(6) 



The pin map fixes the system into a single pure state, 
which means that the state of the environment is fixed 
into a single state as well; the pin map decorrelates the 
system from the environment. Once the pin map, 0, is 
applied, the system is prepared in the various different 
input states by applying local 
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p(m) 



(7) 



As seen the last equations, the advantage of stochastic 
preparation method, as described here, is that the state 
of the environment is a constant for any input of the 
system. This constancy of the state of the environment 



^ Generally we dot not require the preparation map to be trace 
preserving. However, the preparation procedure used here, the 
stochastic procedure, is trace preserving, hence we need to worry 
about this issue. See [8] for details. 



is necessary to characterize the dynamical mechanisms 
properly [7]. 

It may seem that stochastic preparation procedure al- 
leviates the experiment from having an environment that 
depends on the preparation procedure. This is the main 
result of this paper, as we analyze this matter in much 
great detail below to show how inconsistencies in prepara- 
tion procedure can lead to strange experimental results. 



C. Quantum process tomography experiment in 

steps 

The basic steps in a quantum process tomography ex- 
periment are broken down below: 

• Just before the experiment begins, the system and 
environment is in an unknown state p^^ . The sys- 
tem and the environment in general are correlated; 
we are no longer making the weak coupling assump- 
tion. 

• The system is altered to a known input state by 
a preparation procedure. The system and environ- 
ment state after preparation is therefore given by 

The input state is given by taking trace with re- 
spect to environment 



p(™) = Tr£ P^") ® 



• The system is then sent through a quantum pro- 
cess. We consider the evolution to be a global uni- 
tary transformation in the space of the system and 
the environment: 



Jjp(m) 



• Finally the output state is observed. Mathemati- 
cally it is the trace with respect to the environment 



)(™) 



Tr. 



If pi™) 



'P 



(8) 



we will call this the stochastic process equation. 

• Finally using the input and the output states, we 
construct a map describing the process is con- 
structed. 

The procedure above is identical to the procedure to find 
a dynamical map, except for the preparation procedure. 
In the next two sections we will analyze quantum process 
tomography with the two preparation procedures dis- 
cussed in last section. The differences between what we 
find here and dynamical maps will be due to the prepa- 



ration procedures (see Sec. IV for a discussion 
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II. QUANTUM PROCESS TOMOGRAPHY 
WITH STOCHASTIC PREPARATIONS 

In this paper we analyze quantum process tomography 
procedures when the stochastic preparation procedure is 
used to generate the input states. The stochastic prepa- 
ration map's action on a bipartite state is given by Eq. 
[7] While the output states corresponding to the input 
is give by the stochastic process equation, Eq. [8] All 
of our examples below will be based on these two equa- 
tions. Furthermore, all of our examples will have the 
same starting point, namely the unknown state before 
preparation 
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' I + C230'2 ® 0-3) . 



(9) 



This is a correlated state, but it is not an entangled state. 
All examples will also have the same global dynamics in 
common, given by 



U = e 



where 



-iHt 



= Y[ {cos (ujt) I (8) I - i sin (cut) aj (g) aj} ,(10) 



(11) 



Let us now delve into constructing our examples. Start- 
ing with an ideal second quantum process tomography 
experiment. 



2a>X 



Figure 1: The eigenvalues of the process map in Eq. |18| are 
plotted as function of 2uLit. As expected the eigenvalues are 
always positive. The stochastic preparation procedure allows 
the experimentalist to prepare any pure state for the system. 
By convexity, then all possible states of the system are in the 
compatibility domain of the process map in Eq. |18| 



Where the inputs states are 

P^''-^^\{l-a,}, pa-+) = l{I + ai}, 

p(2,+)= + p(3. + )= l{]I + ^3}. (15) 

Now each input state is sent through the quantum pro- 
cess. We will use the following unitary operator to de- 
velop the system and environment. Noting the duals for 
the input states in Eq. [15] are 





p(i,-) 


-\^^ 


- (Ti 


- 0-2 - (T3}, 


A. Ideal stochastic preparation 




-> 








p(i.+) 




- 0-2 - 0-3}, 


Let the pin map, Q be 


p(2, + ) 




p(3- 





e 



(12) 



where \4>) is a pure state of the system. The prepara- 
tion of p^^ in Eq. [9] with this pin map leads to 



(16) 



We can calculate the output states Eq. [8] The output 
states are 



(13) 



yielding the initial state 10) (01 for the system qubit and 
a completely mixed state for the environment qubit. The 
next step is to create the rest of the input states using 
maps fi'^™''. In this case, the fixed state \4>) (0| can be 
locally rotated to get the desired input state (given 
in Eq. 







- cos^(2wi)a-i}. 






+ cos^(2a;<)(Ti}, 






+ cos^(2u;t)(T2}, 






-1- cos^(2u;t)(T3}. 



15) 



10) (0 



(17) 



The linear process map is constructed using Eq. |3] the 
duals in Eq. 16 and the output states: 



(14) 



where m = {(1, -), (1, -h), (2, -h), (3, -h)} and V^"*) are 
local unitary operators acting on the space of the system. 



A. = 



1 + C2 2C2 

1-C2 

1-C2 

2C2 l-l-C^ 



(18) 
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where C = cos{2u>t). The eigenvalues of the process map 
are plotted in Fig. [l] The map above is completely pos- 
itive and can be tested for linearity. This would be the 
result of an ideal quantum process tomography experi- 
ment (see Sec. Ill A for such an experiment). The key 



thing note above is that the state of the environment is 
completely constant throughout the problem. Below we 
will show that this need not be the case all of the times. 



B. Negative maps due to control errors 



The process map turns out to be rather complicated, 
and its eigenvalues are even more complicated looking. 
Therefore, we do not write them down, instead we have 
plotted the eigenvalues as function of iujt. We take the 
value for the error to be e = 0.1 for the plot in Fig[2j One 
of the eigenvalue in Fig. |2]is negative for certain times. 
This shows yet another cause for negative eigenvalues in a 
process map. The negative eigenvalues here have nothing 
to do with the initial correlations between the system and 
the environment. The negative eigenvalues are attributed 
to poor control in the preparation procedure. 



Suppose the initial state is prepared well using a pin 
map 



(19) 



After obtaining this state, the other input states are pre- 
pared by local rotations. Let us consider the case where 
one of the rotation is not perfect. 



1 
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{VT^e\l)-VT+~e\0)), (20) 



where e is taken to be a small positive real number. We 
introduced a small error for the preparation of P^'^'~\ 
but we have kept the error simple by not giving it an 
additional phase, i.e. keeping e to be real. For simplicity 
we have also assumed that the error of this sort occurs 
in the preparation of only one state. 

Let us now pretend that we are not aware of this error. 
Then in reality we have the following set of input states 







f £(73 - - e^CTi}, 












p(2, + ) 


= -{t 

2^ 






p(3, + ) 


2^ 




(21) 



Let us use the same unitary evolution as before given 
in Eq. [TO] The output states corresponding to the last 
three i npu t states in Eq. [21] are the same as before, given 
by Eq. [l7j For the input state P^^'~\ the corresponding 
output state is follows 





^^<- 


h (e(73 - Vl - eVi) cos^{2ujt)}, 






h cos^{2u;t)ai}, 






h cos^{2u;t)a2}, 






h COS^(2cjt)cr3}. 



(22) 



Using these output states, the duals given by Eq. [16] and 
Eq. [3| we can find the process map. 




lojt 



Figure 2: The eigenvalues of the process map found using 
the output states in Eq. [22] and the duals in Eq. |16| One 
of the eigenvalue is negative for certain times; we have taken 
e — 0.1. The negativity is due to the errors in the unitary 
operation implemented to prepare one of the input states. 



Now we have discussed several scenarios that can lead 
to negative eigenvalues for a process map. Two com- 
ments are in order at this point. When the prepared 
input states are not pure and the process map has nega- 
tive eigenvalues, then one should be weary of initial cor- 
relations with the environment. To further check this, 
the process map should be tested for linearity by send- 
ing several additional input states through the process. If 
the process map predicts the output states properly, then 
one can be confident that the system is initially corre- 
lated with the environment. If the process map does not 
predict the output states correctly (non- linear behavior) , 
then there are additional problems with the experiment, 
including the possibility of poor control in the prepara- 
tion procedure. In the case where the input states are 
pure and the process map has negative eigenvalues, the 
negativity can only come from either inconsistencies in 
the preparation procedure or poor preparation control. 



C. Mixed input states 

Demanding pure input states in a quantum process to- 
mography experiment guarantees initially uncorrelated 
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state of the system. Though if it is not possible to pre- 
pare pure states, we can still determine the process. Be- 
low we discuss what happens when the input is not a 
pure state. We will show that this creates an uncertainty 
about whether the system is correlated with the environ- 
ment or not. Yet we show that regardless of this uncer- 
tainty, the quantum process tomography experiment can 
be carried out consistently. 



1. Uncorrelated 

Suppose the initial state is prepared by a pin map lead- 
ing to 

e(/^) = (pP(^'+) + (l-p)p(3-)^ (23) 
= ^{I + pas}®!, (24) 

where << p < 1. The rest of the input states are 
prepared by rotations. Then the input states are 

p(i-) = i{I-pai}, P(i'+' = ^{l + pai}, (25) 
P(2,+)^ + p(3.+)^i{i + p^3}. (26) 

For the unitary operator in Eq. |10[ the corresponding 
output states will be 



(1,-) 




- pcos^{2ujt)ai}, 








\- pcos^{2ujt)(Ti}, 




(2, + ) 








(3, + ) 




|-pcos^(2tjt)cr3}. 


(27) 



The only change that we have to make to find the 
process map is to define a dual proper set, in this case 



1 



P^i--^ = ^{pl-'Ji-a2-a3}, 
2p 

P(i'+) = ^{pl + a^-a^-as}, 
2p 

P(2'+) = ic72, P(3'+) = -a,, 
p p 

We can find the process map using Eq. |3] 



(28) 



/l + C^ 2C^ 

. I I , 

A™.i - - Q Q ^ _ ^2 I ' 

\ 2C^ 1 + C2 

where C = cos(2a;i). The process map here case turns 
out to be the same as in Eq. [l8] This is to expected; the 
process map for stochastic preparation is defined over set 
of all states, including the input states above. 



2. Correlated 



The downside of course is that there is no way to dis- 



tinguish the states in Eq. 23 from the following state 



= ^{I'»I + PCr3 + C23Cr2 0Cr3}. (30) 



Unlike in Eq. 25 this is a correlated state. Yet, both 
total states have the same reduced state for the system 
part. 

Even in this case we can find the process map properly 
if we use the correct dual set given in Eq. [28j The process 
map in that case will be 



/ 1 + C2 -C23CS' 2C2 \ 
-C23CS' 



-C23CS 
V 2C^ -C23CS 





1 + 



where C = cos(2wt) and S — sin(2aji). Which the same 
as the dynamical map calculated in (Eq. 6) [H]. In both 
of the examples above, had we assumed that the input 
states were close enough to the pure states we desired, 
and used the dual set give by Eq. [I6] then the process 
map would contain an error. 




Figure 3: The eigenvalues of the dynamical map in Eq |31| are 
plotted as function of 2u!t. One of the eigenvalue is negative 
for certain values of cut; we have taken C23 — 0.5. The negative 
eigenvalue is due to the initial correlations between the system 
and the environment. 



These two examples illustrate the contrary argument 
to what some have suggested, that only the completely 
positive process maps should be considered physically 
valid. We showed in the second example that, one can 
obtain a not-completely positive process map in a con- 
sistent fashion. The real issue is that the correlations 
with environment are not convenient for experimental 
purposes. But that does not mean that one should fix 
a not completely positive process map to a completely 
positive process map with numerical methods |35l 136] . 
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D. Multiple stochastic preparations 

The pin map used in stochastic preparation must be 
used consistently. Let us show an example of what hap- 
pens when two stochastic preparation procedures are 
used. Consider a quantum process tomography experi- 
ment where the following linearly independent states are 
used to span the space of the system 



pi 



p(l+)^ p{3+)_ 



(31) 



These states form the linearly independent set^ {I, CTj}; 
which is different from the set used in previous examples 

Consider the following two qubit state as the available 
state to the experimenter e± t — 0-: This is the same 
state we used in Eq. [9] for the previous example. 

Let the pin map, G, for our example be 



e 



(32) 



where |(/)) {(f>\ is a pure state of the system. The prepara- 
tion of p^^ with this pin map leads to 



(33) 



yielding the initial state \(p) {<j)\ for the system qubit and 
a completely mixed state for the environment qubit. The 
next step is to create the rest of the input states using 
maps ri'™'. In this case, the fixed state 10) {(f>\ can be 
locally rotated to get the desired input state (given 
in Eq. 



15) 



(34) 



where m = {{l,-),{l,+),{2,+),{3,+)} and V^"''> are 
local unitary operators acting on the space of the system. 

Once again, we take the same initial state and unitary 
operators as in the last example. Now suppose the pin 
map given in Eq. [32] is used to prepare the state 10) (4>\, 
and then by local transformations p(^'+\ p(2,,+)^ and 



p(3,+) are prepared. Finally the mixed state is prepared 
by letting p(3+) decohere. This is a different pin map 
than the one in Eq. [32] so we are using two pin maps to 
prepare two sets of input states. The unitary operator, 
we are using, is often called the swap gate, because it 
swaps the states of two qubits with the period of i = ^ . 



Then at t ■ 



-j^, the total state will be 



p-- 1^^) = i{I (8) I + Qjl (g) CTj}. 

The state of the system has fully decohered. 

The corresponding output states for the input states 
above are found using Eq. [8] The state of environment 
in that equation for input = jl is p^ — ^{I + ajaj}, 
while for the other inputs the state of the environment 



(35) 



^I. Then the corresponding output states are 

Q^'^ = ^{I + sin2(2u;t)a3}, 
Q(''+) = ^{I + cos2(2wt)ai}, 



-{I + cos'^{2ujt)a2}, 



Q^''+) = ^{I + cos2(2wt)a3}. 



(36) 



The last three output states are the same as in the last 
example, but the fourth one is different. 

Now suppose we calculate the output state correspond- 
ing to the input P^~^\ By linearity we have P^^^) = 
I — p(+^^. If the process is linear, then the output state 
for this input state is given by 



Q 



(-1) - 



2Qm - Q(+i) 

i{I + 2sin2(2w<) cTg - cos2(2wt)cri}. (37) 



This state is not physical for certain times, and therefore 
the linearity and the positivity of the process is violated. 

The process has not changed from the last example, 
only the method of determining the process has. Thus 
when the stochastic map is not used consistently, the 
process map can behave nonlinear. For completeness we 
find the process map for this example, 




' C'^ + 25-2 





2C2 


{l-i)S^ 1 + C^-2S^ 



(38) 



r 



where C — cos{2Lut) and S — sm{2ujt). The eigenvalues 
of this process map are negative, as seen in Fig. ]4] The 
fact that the state of the environment is not a constant 
for all of the input states leads to the negativity and 



non-linearity of the process map. 

The origin for the negative eigenvalues here is com- 
pletely different than the negative eigenvalues is not due 
to the initial correlations between the system and the en- 
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2wt 



The experimentalist find the resulting process map to 
be completely positive without any numerical corrections. 
This in principle is an ideal quantum process tomography 
experiment. The absence of negativity from the process 
map suggests that the system of interest was not in any- 
way correlated with the surrounding degrees of freedom. 

However to be sure, one should check this process 
map for linearity by preparing arbitrary initial states and 
putting them through the swap procedure. If the exper- 
imental result matches the theoretical predictions of the 
process map, then the results of this experiment fit the 
analysis of the constructed example in Sec. II A 



Figure 4: The eigenvalues of the process map in Eq. [38] plot- 
ted as function of 2Ljt. One of the eigenvalue is negative for 
certain values of cut. The negativity is due to the inconsis- 
tency in the preparation procedure, and not due to the initial 
correlations with the environment. 



vironmcnt. Since none of the inputs were correlated with 
the environment here, negativity arises from the incon- 
sistencies in the preparation procedure. 



While this example may not seem realistic, the point 
regarding inconsistency arising from multiple pin maps 
still stands. In realistic cases, the trouble may not be 
seen so easily, due to complicated interactions with the 
environment. We will discuss an experiment where mul- 
tiple stochastic preparation procedures are implemented 
in Sec. HITdI 



III. ANALYSIS OF EXPERIMENTS 

In this section we analyze four quantum process to- 
mography experiments that are analogues of the four con- 
structed examples of the last section. Each experiment 
is a well performed experiment in its own right. We an- 
alyze each experimental procedure and comment on why 
the result was not completely positive (when applicable) . 



A. Quantum process tomography of 
Josephson-phase qubit and two-level state 

In this experiment |27j , a Josephson-phase qubit along 
with a more stable but less controllable two-level state 
acting as memory are examined. The Josephson-phase 
is initially prepared into the ground state with high fi- 
delity. Which is then prepared into various excited states 
via unitary transformations. Finally it is 'swapped' with 
the memory qubit and then some time later swapped 
back. Quantum process tomography of this procedure is 
implemented. By doing so, they observe that the deco- 
herence of the Josephson-phase qubit can be hindered by 
transferring the information to the stable two-level state. 



B. Quantum process tomography solid-state qubit 

In this experiment [28 , a solid-state qubit is put 
through various gates. Quantum process tomography of 
these gates is performed. They perform quantum pro- 
cess tomography on three operations; identity, rotations 
along a;— plane, and rotations along y— plane. The results 
are excellent, with very little errors. 

For the case where identity operation is made, the pro- 
cess map is nearly identity as expected. This means that 
the input states are well prepared, and are nearly pure. 
For the other transformations, the errors are attributed 
to errors in preparation procedures and subsequent mea- 
surements. 

The process map obtained in each case, however has 
negative eigenvalues. Since the negativity is along the 
same magnitude for the case when no operation is made 
and when rotations along x and y— planes are imple- 
mented, we can presume that the errors must be due 
to the errors in unitary operations used to prepare input 
states. This is along the same lines as the example in 
Sec. HTbI 



C. Quantum process tomography in NMR 

In this experiment f25], the system that is studied is 
an electron configuration formed in a nitrogen vacancy 
defect in a diamond lattice. The quantum state of the 
system is given by a spin triplet {S = 1). Again we will 
write the initial state of the system and environment as 

The system is prepared by optical pumping, which re- 
sults in a strong spin polarization. The state of the sys- 
tem is said to have a 70% chance of being in a pure state 
10). Or more mathematically, the probability of obtain- 
ing isTr[|0) (0|p^^] = O.7. 

Since the population probability is high, an assump- 
tion was made that the state of the system can be sim- 
ply approximated as a pure state {(j)) ((/)|. From this initial 
state, different input states can be prepared by suitably 
applying microwave pulses resonant with the transition 
levels. After preparation, the system is allowed to evolve, 
and the output states are determined by quantum state 
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tomography. With the knowledge of the input state and 
the measured output states, a linear process map is con- 
structed. 

It was found that the linear process map has negative 
eigenvalues, so the map was "corrected" using a least 
squares fit between the experimentally determined map 
and a theoretical map based on Hermitian parameteriza- 
tion |35j . while enforcing complete positivity. 

However, if we do not regard the negative eigenvalues 
of the map as aberrations, then we should consider the 
assumptions about the preparation of the system more 
carefully. The assumption about the initial state of the 
system is: 



(39) 



This is in effect a pin map. Along with the pin map, 
the stochastic transformations are applied on the initial 
state to prepare the various input states; this is identical 
to the stochastic preparation method discussed in Sec. 

mn 

It is clear that the pure initial state assumption is un- 
reasonable given our knowledge now of how the process 
is sensitive to the initial correlations between the system 
and the environment. In effect the action of the pin map 
in this experiment is not perfect, and the pin map can 
be ignored. Then the process equation is: 



Q 



(m) 



TrB[C/r!("V'^^C^1 



(40) 



where fi^'") is the stochastic transformation that prepares 
the mth input state. In this experiment, fi^™) is noth- 
ing more than a unitary transformation satisfying 
^(™) 10) = |^/;(™)), where \^p^"^^) is the desired pure mth 
input state. 

We can write the unitary transformation for a two-level 
system as: 

y('") = |V-(™))(0| + |^M)^^^| (41) 

where (V'^") IVA"^) = (<^I<?^-l) = 0. This defines as a 
transformation from the basis {10)} to the basis {IV'!"'')}- 
The equation for the process becomes: 



)(™) — 



+Trs [c/|^i'"^)((/)i|p'5^|<^) (V-f^'lC/t 

+Tr£ [c/|^('"))(0|p^^|0^)(vA"^|C/^ 

+Tr, [c/|^i"))(0^|p^^|0^)(^i'"V^ 

Therefore, since {(plp^^ = 0.7 to first approxima- 
tion, the process is a linear mapping on the states 
|^(™)) ('0(™)|. However, it is clear that if all terms are 
included, the process is not truly linear in the states 
|^(m)) (^(m)|^ The negative eigenvalues are therefore a 
result of fitting results into a linear map when the pro- 
cess is not truly represented by a linear map. 



There is no way to tell if the system in the experiment 
above was initially correlated with the environment or 
not, as shown in Sec. |II C[ However the negativity in the 
process map obtained that the system may have been 
correlated initially. If this was the case the obtained map 
may be linear and well behaved. Fixing it numerically to 
make it completely positive may in reality make it behave 
in non linear fashion. 



D. Quantum process tomography of motional 
states 

In this experiment [26], quantum process tomography 
of the motional states of trapped ^^Rb atoms in the po- 
tential wells of a one dimensional optical lattice is per- 
formed. Only two bound bands are considered, which 
are labeled as states |0) and |1). The states are prepared 
stochastically. An initial state of the system is the ground 
state \g) {g\, and from it states |r) (r|, \i) and the fully 
mixed state |l are prepared. The states |r) and \i) stand 
for the real and imaginary coherence states, which are 
prepared by applying appropriate unitary transforma- 
tions on the ground state. This is achieved by displacing 
the lattice for the real coherence state and for the imag- 
inary coherence state, a quarter-period delay is added 
after the displacement. The identity state is prepared by 
letting a superposition state decohere. The input states 
are allowed to evolve and the output states are determine 
by quantum state tomography. 

The process map is found following the usual proce- 
dure laid out in Sec. |Tj Since there are particles lost to 
the neighboring cells, the map is not required to be trace 
preserving. Based on this loss they also argue that the 
map can pick non-physical behavior (not-completely pos- 
itive). The map is forced to be "physical" (completely 
positive) by using the maximum likelihood method |37l . 

In our terminology, the states prepared are 
p(2,+)^ p(3,-)^ ^j-^j ij r^Yie three projective states are 
prepared by a single consistent stochastic preparation, 
while the fully mixed state is prepared by letting the 
state p(i'+) decohere. Which means an additional pin 
map is applied to prepare one of the input states. As 
we saw in the example in Sec. |IID| this can lead to a 
not-completely positive and non- linear process map. 

Furthermore, the input states have varying values for 
polarization. Their data is listed in the table below. 



Pg Pi Pr Pi 



0.90 0.60 0.69 0.69 
P(3.+) 0.10 0.40 0.31 0.31 
P(i'+) 0.82 0.59 0.85 0.63 

0.84 0.58 0.64 0.37 



where pj are the experimentally prepared states pro- 
jected onto projectors P^"^\ The polarization of the 
imaginary state in the cr2 direction is very low. This could 
mean that it is correlated with the environment, while 
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the polarization of the ground state along the negative 
(73 direction is almost unity, meaning it is only weakly 
correlated at best. 

Here, we have two potential causes for the process map 
to have negative eigenvalues. The first problem is with 
the experimental procedure; applying multiple stochas- 
tic preparations that may affect the state of the environ- 
ment differently. The second problem may be unavoid- 
able; since it is extremely difficult to prepare pure states 
in a setup like this. Even if the negative eigenvalues are 
due to the initial correlations, we do not have a prescrip- 
tion to obtain a process map in a consistent fashion. For 
the last example shown in Sec. [IT] we assumed that the 
initial correlations were constant throughout the process. 
This does not seem to be the case here. Since the ground 
state is almost pure, while the imaginary state is clearly 
not pure. Therefore the correlation with environment for 
these two inputs must be different as in the case of the 
example in Sec. |fID[ 

IV. CONCLUSION 

The examples with multiple stochastic preparations (in 
Sec. |ff D| and projective preparations (in Sec. V.A) 
in [5]) look very similar. For multiple stochastic prepa- 
rations, the state of the environment depends on the 
stochastic map. If we think of each projective prepara- 
tions as an independent stochastic preparation then the 
state of the environment, in the example (in Sec. V.A) 
in depends on the preparation map. Then are the 
two situations the same? Let us show below that the two 
situations are fundamentally different. 

In the weak coupling limit, the projective preparation 
procedure does not play an important role. As we saw in 
the last section, when the system and the environment 
are initially uncorrelated, the projective preparation pro- 
cedure will have no affect on the state of the environment. 
This is because the state of the environment is affected 
only indirectly due to the initial correlations between the 
system and the environment. Thus the projective prepa- 
rations, in the weak coupling limit, yield a linear process 
map. 

This is not the case in the multiple stochastic prepara- 
tion example. When multiple stochastic maps are used 
to prepare different input states, the inconsistencies do 
not stem from the initial correlations between the sys- 
tem and the environment. In fact, in our example all 



input states are initially uncorrelated from the environ- 
ment. The inconsistencies arise from the preparation pro- 
cedures themselves, leading to different states of the envi- 
ronment for different input states. These inconsistencies 
will be absent in the case where the system develops in a 
closed form, since the system does not feel the presence 
of the environment during the quantum process. Hence 
in the "weak interaction limit" multiple stochastic prepa- 
rations yield a linear process map. 

Lastly, note that each of the process map found in this 
section, given by Eqs. 38 and (in Sec. V.A) in [5] are 
different from each other. This clearly shows that the 
preparation procedures play a non-trivial role in open 
quantum system experiments. The preparation proce- 
dure is the only thing that distinguished each case. For 
the case where no preparation procedure is applied, i.e. 
the case of the dynamical map (Eq. 6 in [9 ), the sit- 
uation is still different. The dynamical map, which has 
negative eigenvalues, is linear and has a valid interpreta- 
tion within the compatibility domain. The process maps 
in Eqs. 38 and (in Sec. V.A) in fSJ do not have any 



consistent interpretation. 

When negative eigenvalues are found in a process map, 
they hint to some problem in the preparation procedure. 
Although this is a bit premature to state at this point; we 
analyze the negative eigenvalues in process map in more 
detail in the next section. 

Our study of quantum process tomography started by 
noting that the dynamical map acting on a system can 
have negative eigenvalues. The dynamical map has neg- 
ative eigenvalues when the system is initially correlated 
with the environment. In the course of our studies of 
quantum process tomography, we showed that the prepa- 
ration procedure cannot be neglected for any quantum 
system that interacts with an environment. These are 
the two major themes discussed in this paper. Though, 
along the way, we presented a method of quantum pro- 
cess tomography that is independent of the preparation 
procedure. The map arising from this procedure lead 
us to an expression that quantifies the memory effect on 
the dynamics of the system due to the initial correlations 
with the environment. Determining the memory effect is 
an important task in coherence control. 
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